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Abstract 

We briefly summarize theoretical methods for carrying out 
QCD calculations to next-to-leading order in perturbation 
theory. In particular, we describe a new general algorithm 
that can be used for computing arbitrary jet cross sections 
in arbitrary processes and can be straightforwardly imple- 
mented in general-purpose Monte Carlo programs. 
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1 Motivations 



During the last fifteen years many efforts have been devoted to carry 
out accurate QCD calculations to higher perturbative orders. These 
calculations are motivated by three main reasons. 

First of all, the comparison between perturbative calculations and 
experimental data allow one to perform precision tests of QCD in 
the strong-interaction processes that involve a large transferred mo- 
mentum Q ll], §]. These tests are essential for measuring the strong 
coupling as(Q) and its running 0] as predicted by asymptotic free- 
dom. Perturbative QCD studies are also important to evaluate the 
background for new physics signals. An outstanding example of that 
is the current investigation [|^ of the discrepancy between the single- 
inclusive jet distribution at large pt, as measured by CDF Q, and the 
QCD predictions. More recently, a renewed interest in perturbative 
calculations has been motivated by phenomenological and theoretical 
models of non-perturbative phenomena (see |^ and references therein). 
Using these models and having under control the perturbative com- 
ponent, one can use experimental data on high-energy cross sections 
to extract information on the underlying non-perturbative dynamics. 

To these aims, calculations at the leading order (LO) of the per- 
turbative expansion in the QCD coupling a;s(Q) are insufficient. In 
fact, just because of its perturbative nature, the running of the QCD 
coupling can be hidden in higher-order corrections. Thus at LO the 
value of as is essentially undetermined and a LO calculation predicts 
only the order of magnitude of a given cross section and the rough 
features of a certain observable. The accuracy of the perturbative 
QCD expansion is instead controlled by the size of the higher-order 
contributions. Any definite perturbative QCD prediction requires (at 
least) a next-to-leading order (NLO) calculation. 

In general, NLO calculations are highly non-trivial. The first bot- 
tleneck one encounters in producing new NLO calculations for a cer- 
tain process is the evaluation of the relevant matrix elements. In 
recent years new techniques [|] have been developed to compute QCD 
Feynman diagrams and most of the one-loop five-point amplitudes 
are now available [0, |^. However, even when the process-dependent 
matrix elements are known, there are practical difficulties in setting 
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up a straightforward calculational procedure. The physical origin of 
these difficulties is in the necessity of factorizing the long- and short- 
distance components of the scattering processes and is reflected in the 
perturbative expansion by the presence of divergences. QCD theorems 
guarantee that these divergences eventually cancel in the evaluation 
of physical cross sections but do not prevent their appearance in in- 
termediate steps. Since single intermediate expressions are usually 
divergent, the numerical implementation of NLO calculations forms a 
second bottleneck. 

The main issue one has to face is thus the following. On one 
side many different NLO calculations (i.e. calculations for different 
observables) for a certain process and, possibly, for many processes 
are warranted. On the other side each calculation is very complicated 
(see also Sect. 

In particular, it is very important to reduce the second bottleneck 
by setting up efficient and simple methods for computing arbitrary 
quantities in a single process. It is even more important to have at 
our disposal simple algorithms for computing arbitrary quantities in 
arbitrary processes. The goal is a universal algorithm that, in princi- 
ple, can be used to construct a general-purpose Monte Carlo program 
(not a Monte Carlo event generator) for carrying out NLO QCD cal- 
culations. Conceptually, such an algorithm could be used in the same 
manner as some universal Monte Carlo event generators (e.g. HER- 
WIG 0): any time one wants to compute a new quantity or to vary 
the experimental cuts, one simply modifies the 'user routine' accord- 
ingly; any time one wants to study a different process, one simply 
enters the corresponding matrix elements. 

A new general algorithm of this type was recently presented . 
It is based on two key ingredients: the subtraction method for the 
numerical cancellation of the divergences among different contribu- 
tions; and the dipole factorization theorems for the universal (process- 
independent) analytical treatment of individual divergent terms. 

In this contribution, after a brief summary of general methods, 
we describe these two ingredients and show some numerical results 
for the specific cases of jets in e~^e~ annihilation and deep-inelastic 
lepton-hadron scattering (DIS). 
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2 NLO QCD calculations 



The general structure of a QCD cross section in NLO is the following 

Here the LO cross section is obtained by integrating the fully 
exclusive cross section da^ in the Born approximation over the phase 
space for the corresponding jet quantity. Let us suppose that this LO 
calculation involves m partons with momenta pk {k = 1, ...,m) in the 
final state. Thus, we write 

a^o = f da^ ^ (2) 

J m 

where the Born-level cross section is: 

da^ = ci$(")({pj) \MU{Pk})\' Ffr\{p,}) , (3) 

and d^^"^^ and A4m respectively denote the full phase space and the 
tree-level QCD matrix element to produce m final-state partons. These 
are the factors that depend on the process. 

The function Fj^"^ defines the physical quantity that we want to 
compute, possibly including the experimental cuts. Note that this 
quantity has to be a jet observable, that is, it has to be infrared and 
collinear safe: its actual value has to be independent of the number 
of soft and collinear particles in the final state. Thus, we should have 
(we refer to for a more detailed formal definition) 

in any case where the m + 1-parton configuration on the left-hand side 
is obtained from the m-parton configuration on the right-hand side 
by adding a soft parton or replacing a parton with a pair of collinear 
partons carrying the same total momentum. 

Efficient techniques, based on helicity amplitudes [|1T| and colour 



subamplitude decomposition [|^, are available for calculating tree- 



level matrix elements. Thus the evaluation of the LO cross section does 
not present any particular difficulty. Even if a^^ cannot be computed 
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analytically (because Aim is too cumbersome or the phase-space cuts 
in Fj"^^ are very involved), one can straightforwardly use numerical 
integration techniques, for instance, a Monte Carlo program where 
the function Fj"^^ is given as 'user routine'. 

At NLO one has to consider the exclusive cross section da^ with 
m + 1 partons in the final state and the one-loop correction da^ to 
the process with m partons in the final state: 

= / rfa^ + / da"" . (5) 

The exclusive cross sections da^ and da^ have the same structure as 
the Born-level cross section in Eq. (0), apart from the replacements 

IMml"^ IMm+il"^ and iTWmP ^ \Mm\fi-ioop)- Here \Mm\fi-ioop) de- 
notes the QCD amplitude to produce m final-state partons evaluated 
in the one-loop approximation. 

The calculation of the loop integral in iMmlfi-ioap) leads to ultra- 
violet, soft and collinear singularities. The ultraviolet singularities can 
be handled in a simple way within the loop corrections by carrying out 
the renormalization procedure. Thus we can assume that the virtual 
cross section in Eq. is given in terms of the renormalized matrix 
element and the ultraviolet divergences have been removed. 

Soft and collinear singularities instead lead to the main problem. 
These singularities do not cancel within the sole da^ and are accom- 
panied by analogous singularities arising from the integration of the 
real cross section da^. In the case of jet quantities, adding the real 
and virtual contribution, these singularities cancel and the physical 
NLO cross section in Eq. is finite. This cancellation is guaranteed 
by the property in Eq. (|[). However, the cancellation mechanism is 
not trivial because it does not take place at the integrand level. 

The two integrals on the right-hand side of Eq. are separately 
divergent so that, before any numerical calculation can be attempted, 
the separate pieces have to be regularized. The most widely used 
regularization procedure (actually, the only regularization procedure 
that is gauge invariant and Lorentz invariant to any order of the QCD 
perturbative expansion) is obtained by means of analytic continuation 
in a number of space-time dimensions d = A — 2e different from four. 
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Using dimensional regularization, the divergences (arising out of the 
integration) are replaced by double (soft and collinear) poles and 
single (soft or collinear) poles 1/e. Thus the real and virtual contribu- 
tions should be calculated independently, yielding equal-and-opposite 
poles in e. These poles have to be combined and, after having achieved 
their cancellation, the limit e — > can be safely carried out. 

In principle this computation procedure does not pose any prob- 
lems. In practice, it is not so. On one side, analytic calculations 
are impossible for all but the simplest quantities because of the in- 
volved kinematics for multi-parton configurations and of the compli- 
cated phase-space cuts relative to the definition of the jet observable. 
On the other side, the use of numerical methods is far from trivial be- 
cause real and virtual contributions have to be integrated separately 
over different phase-space regions and because of the analytic contin- 
uation in the arbitrary number d of space-time dimensions. 

The most efficient solution to this practical problem consists in 
using a hybrid analytical/numerical procedure: one must somehow 
simplify and extract the singular parts of the cross section and treat 
them analytically; the remainder is treated numerically, independently 
of the full complications of the jet quantity and of the process. 



2.1 General methods and algorithms 

There are, broadly speaking, two general methods for doing that: the 
phase-space slicing method and the subtraction method. Both the 



slicing |I3| and the subtraction [14| methods were first used in the 



context of NLO calculations of three-jet cross sections in e'^e~ annihi- 
lation. Then they have been applied to other cross sections, adapting 
the method each time to the particular process. Only recently has 
it become clear that both methods are generalizable in a process- 
independent manner. The key observation is that the singular parts 
of the QCD matrix elements for real emission can be singled out in a 
general way by using the factorization properties of soft and collinear 
radiation |T^. Owing to this universality, the two methods have led 
to general algorithms for NLO QCD calculations. 

In the context of the phase-space slicing method, an algorithm has 
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been developed for jet cross sections in lepton and hadron collisions 
IC, 17 . The generalization of this method to include fragmentation 



functions and heavy flavours is considered in Refs. [T^ . 

As for the subtraction method, two approaches are available for 
setting up general algorithms. The 'residue approach' introduced in 
Ref. ^ 



20 1 has been further generalized in Refs. El^, p2 



The dipole 

formalism [24| has been completely worked out in Ref. ||1CI|| . 

The advent of these algorithms has made feasible NLO QCD calcu- 
lations for multi-jet cross sections. Monte Carlo programs have been 
constructed for most of the physical processes that involve four par- 
ticles at LO. For five-particle processes, the three-jet cross section in 
hadron collisions in the simplified case of pure-gluon subprocesses is 
available | ^ 

the full 



, as is the four-jet cross section in electron-positron anni- 
hilation in the approximation of large number of colours [EB 
QCD results are expected to appear soon. 

We refer to Sect. 12.2 of Ref. [l^ for a discussion of the comparison 
among different general methods for NLO calculations. In the rest 
of this contribution we describe the approach, based on the dipole 
formalism. 



3 The subtraction method 

The general idea of the subtraction method is to use the identity 





da^ - da^ 


+ f da^+ 1 


Jm+l 




Jm+l Jm 



which is obtained by subtracting and adding back the same quan- 
tity da^. The cross section contribution da"^ has to fulfil two main 
properties. 

i) Firstly, it must be a proper approximation of da^ such as to 
have the same pointwise singular behaviour (in d dimensions) as da^ 
itself. Thus, da^ acts as a local counterterm for da^ and one can 
safely perform the limit e ^ under the integral sign in the first 
term on the right-hand side of Eq. @. This defines a cross section 
contribution 0-^LO{m+i} ^j^]^ ^ _|_ i-parton kinematics that can be 
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integrated numerically in four dimensions: 



^NLOim+l} 



I Uda"") -(da^) 1. (7) 



a) The second property of da^ is its analytic integrability (in 
d dimensions) over the one-parton subspace leading to the soft and 
collinear divergences. In this case, we can rewrite the last two terms 
on the right-hand side of Eq. as follows 

^NLO{m}^ j L j 
J m \_ Jl 

Performing the analytic integration da^, one obtains e-pole contri- 
butions that can be combined with those in da^ , thus cancelling all the 
divergences. The remainder is finite in the limit e — and thus defines 
the integrand of a cross section contribution (j^-f-oim} ^^^ith m-parton 
kinematics that can be integrated numerically in four dimensions. 
The final structure of the NLO calculation is as follows 

^NLO ^ ^NLO{m+l} _^ ^NLO{m} ^ 

and can be easily implemented in a 'partonic Monte Carlo' program, 
which generates appropriately weighted events with m + 1 and m final- 
state partons. 

Note that, using the subtraction method, no approximation is ac- 
tually performed in the evaluation of the NLO cross section. Rather 
than approximating the cross section, the subtracted contribution da^ 
defines a fake cross section that has the same dynamical singularities 
as the real one and whose kinematics are sufficiently simple to permit 
its analytic integration. 

The real cross section contribution da^ has the following general 
structure 

da« = c^^^-^+i) \Mm+i{{Pk})\' F!r+'\{p,}) , (10) 

where and |A^m+ip depend on the process and i^j™"^^^ de- 

pends on the quantity we want to compute. Obviously, for any given 
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da^ one can try to construct a corresponding da^ by properly approx- 
imating d(^''^^^\ \M.m+i\^ and p^^^^\ It is less obvious that one can 
use the subtraction method to compute arbitrary quantities in a given 
process, because one needs a fake cross section da"^ that depends only 
on the process and, hence, is independent of the actual definition of 
the jet function p^^^\ It is still less obvious that one can use the 
subtraction method to construct a universal algorithm for computing 
arbitrary quantities in arbitrary processes. To this purpose the fake 
cross section da^ also has to be somehow independent of M.m+i- 

Our method to achieve this generality is based on the dipole for- 
malism. 



4 Dipole formalism and universal subtrac- 
tion term 

4.1 Soft and collinear limits 

The starting point of the dipole formalism are the soft and collinear 
factorization theorems for the QCD matrix elements. According to 
these theorems, the singular behaviour in d dimensions of a generic 
tree-level matrix element •••,Pm+i) with m + 1 final-state 

partons can be obtained by means of factorized limiting formulae that, 
respectively in the soft (when the parton momentum pj vanishes) and 
collinear (when the parton momenta pi and pj become parallel) re- 
gions, have the following structure 

)P ^ |-Mm(Pl,---,Pm+l)P ®c J^(Pj) , (11) 

\M 

m+1 

)| ^ \Mm{Pl,---,Pj+Pi,---,Pm+l)\ ®hPij 

(12) 

The notation in Eqs. (|n|,^) is symbolic (see Ref. for more details) 
but sufficient to recall their main features. 

The contributions A^m on the right-hand sides are the tree-level 
matrix elements to produce m partons and are respectively obtained 
from the original m + 1-parton matrix element by removing the soft 
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parton pj or combining the two collinear partons pj and pi into a 
single-parton momentum. 

The other contributions on the right-hand sides are responsible for 
the soft and colhnear divergences. The factor J'^{pj) in Eq. (|ll|) is 
the eikonal current for the emission of the soft gluon pj, and Pij is the 
Altarelh-Parisi sphtting function. These factors are universal: they 
do not depend on the process but only on the momenta and quantum 
numbers of the QCD partons in A4rn- In particular, J'^{pj) depends on 
the colour charges of the partons in A^^, and Pij depends on their he- 
licities. Because of these colour and helicity correlations (symbolically 
denoted by ®c and Eqs. are not real factorized expres- 



sions. Moreover, there is another important reason, due to kinematics. 



why Eqs. ([TT| , p!2|) cannot be regarded as true factorization formulae but 
rather as limiting formulae. Indeed, the tree-level matrix elements in 
Eqs. (|Tl],|T2|) are unambiguously defined only when momentum conser- 
vation is fulfilled exactly. Since, in general, the m + 1- parton phase 
space does not factorize into an m-parton times a single-parton phase 
space, the right-hand sides of these equations are unequivocally de- 
fined only in the strict soft and collinear limits. 

Owing to their universality, the limiting formulae (0,0) can be 
used to approximate the matrix element |A^m+iP in Eq. ( [T0| ) and thus 
to find a fake cross section da^ that matches the real cross section da^ 
in all the singular regions of phase space. However, the implementa- 
tion of Eqs. (11], 12) in the calculation of QCD cross sections requires 
a careful treatment of momentum conservation away from the soft 
and collinear limits. Care also has to be taken to avoid double count- 
ing the soft and collinear divergences in their overlapping region (e.g. 
when Pj is both soft and collinear to pi). The use of the dipole factor- 



ization theorem introduced in Ref. I^J] allows one to overcome these 
difficulties in a straightforward way. 

4.2 Dipole formulae 

The dipole factorization formulae have the following symbolic struc- 
ture 

\Mra+l{Pl,---.Pm+l)? =\M^{pi,...,P^)\^ ®Vij + ... . (13) 
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The dots on the right-hand side stand for contributions that are not 
singular when Pi ■ Pj — >■ 0. The dipole sphtting functions Vij are 
universal (process-independent) singular factors that depend on the 
momenta and quantum numbers of the m partons in the tree-level 
matrix element Colour and helicity correlations are denoted 

by the symbol ®. The set pi,...,pjn of modified momenta on the 
right-hand side of Eq. (|T3|) is defined starting from the original m + 1 
parton momenta in such a way that the m partons in |A^mP are 
physical, that is, they are on-shell and energy-momentum conservation 
is implemented exactly: 

= , Pl + ...+Pm=Pl + ■■■+Pm+l ■ (14) 

The detailed expressions for these parton momenta and for the dipole 
splitting functions are given in Ref. . 



Apart from the presence of colour and helicity correlations, Eq. ([13|) 
can be considered as a true factorization formula because its left-hand 
and right-hand sides live on the same phase-space manifold. Equa- 
tion (|T^ indeed guarantees that exact kinematics are retained in the 
definition of the m-parton configuration {pi, ...,pm}- These m par- 
ton momenta depend on pi and pj in such a way that in the soft and 
coUinear regions the m-parton configuration become indistinguishable 
from the original m + 1-parton configuration. Correspondingly, the 
dipole splitting function Vjj is defined in order to coincide with the 
eikonal current and with the Altarelli-Parisi splitting function respec- 
tively in the soft and collinear limits. 

It follows that Eq. ([13D provides a single formula that approxi- 
mates the real matrix element |A1m+ip for an arbitrary process, in 
all of its singular limits. These limits are approached smoothly, thus 
avoiding double counting of overlapping soft and collinear singular- 
ities. The exact implementation of momentum conservation makes 
possible this smooth transition and the extrapolation of the limiting 



formulae (|Tl],|T2D away from the soft and collinear regions. 
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4.3 Universal subtraction term 



These main features of the dipole formulae allow us to construct a 
universal subtraction term with the following form 

da"" = \Mm{{Pk})\' ® Vi, Ff ^(fe}) . (15) 

Note that the only dependence on the jet observable is in the jet- 
defining function Fj"^^ and the only dependence on the process is in 
the tree-level matrix element |7Vlmp. These are the same m-parton 
functions as enter in the calculation of the Born-level cross section of 
Eq. (Ip. The only other ingredients needed to construct da"^ are the 
dipole splitting functions, which are completely process-independent 



and given once and for all |T^. This specifies the universal character 
of Eq. ([T5|) : the fake cross section da^ used for the NLO calculation 
is straightforwardly obtained in terms of the sole (process-dependent) 
information that is necessary for the corresponding LO calculation. 

Having the subtraction term in the explicit form (^), we can dis- 
cuss how it fulfils the properties i) and ii) listed in Sect. |^. As for 
the property i), we note that there are several dipole terms on the 
right-hand side of Eq. (|1^). Each of them mimics one of the m + 1- 
parton configurations in da^ that are kinematically degenerate with 
a given m-parton state. Any time the m + 1-parton state in da^ ap- 
proaches a soft and/or collinear region, there is a corresponding dipole 
factor in dcr^ that approaches the same region with exactly the same 
probability as in da^. The equality of the two probabilities directly 
follows from (^) and from the limiting behaviour in Eqs. (^,TT,T^) of 
the cross section factors on the right-hand side of Eq. (0). In this 
manner da^ acts as a local counterterm for da^. Note, in particular, 
that the cancellation mechanism is completely independent of the ac- 
tual form of the jet-defining function and works for any jet observable 
(i.e. for any quantity that fulfils Eq. (^). 

As for the property ii) , we start by noting that da^ (likewise da^) 
depends on the m + 1 parton momenta pi, ...,Pm+i- However, having 
introduced the modified momenta pi,...,pm, for each dipole term in 
Eq. ([ISD we can define a one-to-one mapping 

{pi,...,Pm+i}^{pi,:.,Pm,Pi+Pj} ■ (16) 
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The key feature of this mapping is that the m modified momenta can 
be chosen in such a way that they obey exact phase-space factorization 
as follows 

c?<l>(™+i)(pi,...,p^+i) =rf$('")(pi,...,p„) + , (17) 

where dip is a single-particle subspace that, for fixed de- 



pends only on the dipole momenta pi and pj [T^. Owing to the exact 
phase-space factorization and to the fact that the fake cross section in 
Eq. (p!5| ) is proportional to the jet quantity calculated from the mod- 
ified m-parton configuration, the integration of the singular dipole 
contributions can be completely factorized (modulo colour and helic- 
ity correlations) with respect to a term that exactly reproduces the 
Born-level cross section: 

/ da^ = [ d^^-\{p,}) \MU{Pk})\' F!r\{Pk}) 

® E / ^<^({?,})fe Vi, = [ da''® I{{pu}) . (18) 

The last factor on the right-hand side of Eq. (|I8]) is defined by 

and contains all the soft and coUinear singularities that are necessary 
to compensate those in the virtual cross section da^ . Owing to the 
convenient definition of the dipole splitting function Vij, it is possible 
to carry out analytically the integration in Eq. (|1^) over the dipole 
phase space in d dimensions. This leads to an explicit and universal 
expression [0] for the factor /, whose e-poles cancel those in the one- 



loop matrix element. 



5 Final results and numerical implemen- 
tation 

The discussion in the previous section shows that, by using the sub- 
traction method and the dipole formulae, one can extract and treat 
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analytically the singular parts of any NLO cross section in a way that 
is independent of the exact details of the observable and of the pro- 
cess. This leaves a remainder that depends on the full complications 
of the jet quantity, but which is finite so that it can be treated either 
numerically or analytically (whenever possible). 

In general, the use of numerical integration techniques (typically, 
Monte Carlo methods) is certainly more convenient. First of all, the 
numerical approach allows one to calculate any number and any type 
of observable simultaneously by simply histogramming the appropriate 
quantities, rather than having to make a separate analytic calculation 
for each observable. Furthermore, using the numerical approach, it 
is easy to implement different experimental conditions, for example 
detector acceptances and experimental cuts. 

In order to summarize the final results of our algorithm and to 
describe their numerical implementation, we start by recalling how 
the LO cross section in Eq. (|^) is evaluated by using a Monte Carlo 
program. One first generates an m-parton event in the phase-space 
region f^$(™) and gives it the weight Then this weighted event 

is analysed by a user routine according to the actual definition of the 
phase-space function Fj^^ and inserted into a corresponding histogram 
bin. 

Following the decomposition in Eq. (^), the NLO cross section 
is obtained by adding two contributions (which are not necessarily 
positive definite) with m-parton (as in the LO calculation) and m + 1- 
parton kinematics, respectively. Unlike the original real and virtual 
contributions, these two terms are separately finite and can directly 
be integrated in four space-time dimensions. 



5.1 The term with m-parton kinematics 

The first contribution is obtained by inserting Eq. (|18D into Eq. (P) 
and can be written as follows 

J m 
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where the master function J-'m{{Pk}) is exphcitly given by 

y'miiPk}) = { \Mrai{pk})\l-loop) + \MU{Pk})\' ^ Ii{Pk}) } • 

^21) 

The first term in the curly bracket is the one-loop renormalized ma- 
trix element for producing m final-state partons. The second term 
is obtained by combining the tree-level matrix element to produce m 
partons and the universal factor I in Eq. (|19|). These two terms are 
defined in c? = 4— 2e dimensions. Owing to the progress made in recent 
years in the analytical techniques for evaluating one-loop amplitudes 
[^], many of them have been calculated. The explicit expression of the 
universal factor I is provided by our algorithm. Thus, one has to carry 
out the expansion in e-poles of the two terms in the curly bracket, can- 
cel analytically (by trivial addition) the poles and perform the limit 
e — > 0. This simple algebraic manipulation is sufficient to construct 
an effective m-parton weight, the master function JF, that is finite. As 
a result, Eq. (|20| ) can be handled by the Monte Carlo program exactly 
in the same way as the LO cross section. 

Note that the two terms on the right-hand side of Eq. ( PT] ) sepa- 
rately depend on the regularization prescription of the soft and collinear 
divergences, namely dimensional regularization. Since different ver- 
sions of dimensional regularization can be used to compute the one- 
loop matrix element, the second term in the curly bracket has to be 
evaluated accordingly. Alternatively, one can fix the latter and use the 
transition rules derived in Ref. ^\ to relate the one-loop amplitudes 



in different dimensional-regularization schemes. 

The necessity to consistently regularize the separately divergent 
components of the cross section is a common feature of any NLO cal- 
culation, independently of the method that is actually used in the 
computation. Failure in the consistent implementation of the regu- 
larization procedure leads to violation of unitarity and, ultimately, to 
an incorrect (although possibly finite) final result. The dipole formal- 
ism is extremely efficient to guarantee unitarity because all the diver- 
gences are isolated in the right-hand side of Eq. (pTl). As explained 



in Ref. |28], for any regularization prescription that is unambiguously 



defined at the level of one-loop matrix elements, one can compute in 
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a simple and consistent way the universal factor I that provides the 
finite and unitary master function J-'. 

5.2 The term with m + 1-parton kinematics 

The NLO contribution with m+ 1-parton kinematics, which is obtained 
by subtracting the fake cross section in Eq. (|T5|) from the real cross 
section in Eq. (p!OD, has the following explicit expression: 

^NLO{m+l} ^ f ^$(m+l) (22) 
Jm+1 

■ f^\Mm+i{{Pk})\' F!r^'\{p,}) - j: \MU{Pk})\' ® • 

The terms in the curly bracket define an effective matrix element that 
is integrable in four space-time dimensions. It follows that the NLO 
matrix element Aim+i, with m -\- 1 final-state partons, can be di- 
rectly evaluated in c? = 4 dimensions, thus leading to an extreme 
simplification of the Lorentz algebra. Knowing the tree-level matrix 
elements and the dipole splitting functions, the Monte Carlo integra- 
tion of Eq. ( P^ ) is straightforward. One simply generates an m + 1- 
parton configuration and uses it to define an event with positive weight 
+ and several counter-events, each of them with the negative 

weight — I^Vlmp (S> Vjj. Then these event and counter-events are anal- 
ysed by the user routine. The role of the two different jet functions 
p[m+i) p{m) of binning the weighted event and counter- 

events into different bins of the jet observable. Any time that the 
generated m + 1-parton configuration approaches a singular region, 
the event and one counter-event fall into the same bin and the cancel- 
lation of the large positive and negative weights takes place. 

6 Monte Carlo programs 

Generalizing the procedure for constructing NLO Monte Carlo pro- 
grams for arbitrary quantities has several advantages. These are prin- 
cipally due to the reduction in the number and complexity of ingre- 
dients that have to be calculated for each new process, and because 
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the d- dimensional integrals only need be done once and can be eas- 
ily checked independently, rather than being buried inside a specific 
calculation. 

Using the general algorithm described in this contribution, we have 
already constructed two Monte Carlo programs (they can be obtained 
from http : // suryal 1 . cern . ch/ user s/seymour/ nlo/) , EVENT2 and 
DISENT. 

EVENT2 computes three-jet observables in e^e~ annihilation. 
In the case of un-oriented three-jet events, this program is comparable 
and in agreement with the program EVENT |^ , which is based on the 
subtraction procedure of Ref. [jl4| and has been used for most of the 
QCD analyses at LEP and SLC [0, As an example we show the NLO 
coefficients for the thrust and C-parameter distributions in Fig. |l|. We 
find that, in general, the numerical convergence of EVENT2 is similar 
to the program of Ref. except close to the two-jet region in which 
ours becomes progressively better. In the case of oriented events i30[| , 
EVENT2 should be compared with a corresponding program, EERAD 
based on the phase-space slicing method. 

DISENT |10, 31| is a NLO program for 2 + 1-jet observables in DIS. 
The program uses the matrix elements evaluated by the Leiden group 
p^ . In Fig. ^ we show as an example the differential jet rate as a 
function of jet resolution parameter fcut, using the k± jet algorithm [0] 
at HERA energies |Q. We see that the NLO corrections are generally 
small and positive, except at very small fcut (where large logarithmic 
terms, —aglog^ fcut, arise at each higher order). In Fig. ||b, we show 
the variation of the jet rate at a fixed fcut with factorization and 
renormalization scales. The scale dependence is considerably smaller 
at NLO. DISENT can be compared with the Monte Carlo MEPJET 



T3. 



[p5|] that uses the phase-space slicing algorithm of Ref. 

The results of the algorithm based on the dipole formalism have 
also been implemented in a program 
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for the calculation of NLO 
QCD corrections to four-fermion final states in e"'"e~ annihilation. 
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THRUST 




Figure 1: Coefficient of (as/27r)^ for the thrust and C-parameter dis- 
tributions. The dotted histograms show the size of the statistical 
errors. 
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Figure 2: Jet cross sections in ep collisions at HERA energies 
{^/s — 300 GeV). (a) The distribution of resolution parameter fcut 
at which DIS events are resolved into (2 + 1) jets according to the k± 
jet algorithm. Curves are LO (dashed) and NLO (solid) using fac- 
torization and renormalization scales equal to Q^, and the MRSD'_ 
distribution functions. Both curves are normalized to the LO cross 
section, (b) The rate of events with exactly (2 + 1) jets at fcut — 0.25 
with variation of renormalization (solid) and factorization (dashed) 
scales. Normalization is again the LO cross section with fixed factor- 
ization scale. 
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7 Summary and outlook 



The calculation of jet cross sections in perturbative QCD requires the 
integration of multiparton matrix elements over complicated phase- 
space regions that depend on the actual definition of the jet observ- 
ables and on the experimental cuts. In general, these phase-space 
integrations can be carried out only by using numerical methods. Be- 
yond LO, however, numerical techniques cannot straightforwardly be 
applied because real-emission contributions and virtual contributions 
are separately divergent. These divergences have to be first regu- 
larized, then evaluated analytically, combined together and cancelled 
before any numerical calculation can be attempted. 

General methods are now available to overcome all the analytical 
difficulties related to the treatment of soft and coUinear divergences in 
NLO calculations. In this contribution we have mainly described one 
of these general formalisms, which has been used to set up an explicit 
algorithm to compute NLO jet cross sections. 

The algorithm combines the subtraction method and the dipole 
formulae to carry out all the analytical work that is necessary to eval- 
uate and cancel the singularities. The final output of the algorithm 
is given in terms of effective matrix elements that can be automati- 
cally constructed starting from the original (process-dependent) ma- 
trix elements and universal (process-independent) dipole factors. The 
effective matrix elements can be integrated numerically or analytically 
(whenever possible) over the available phase space in four dimensions 
to compute the actual value of the NLO cross section. If the numerical 
approach is chosen, Monte Carlo integration techniques can be easily 
implemented to provide a general-purpose Monte Carlo program for 
carrying out NLO QCD calculations in any given process. 

The simplified discussion of the algorithm presented in this con- 
tribution directly applies to processes, like e+e" — > n jets, in which 
there are neither initial-state hadrons nor identified hadrons in the 
final state. However, the formalism and the algorithm are completely 
general in the sense that they apply to any jet observable in a given 
scattering process as well as to any hard-scattering process. Full de- 
tails and explicit results for lepton-hadron and hadron-hadron colli- 
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sions and for fragmentation processes are given in Ref. |T^. 

At present, next-to- next-to-leading order (NNLO) QCD calcula- 
tions are feasible only for some fully inclusive quantities In these 



cases one considers all possible final states and integrates the QCD 
matrix elements over the whole final-state phase space. Thus one can 
add real and virtual contributions before performing the relevant mo- 
mentum integrations in such a way that only ultraviolet singularities 
appear at the intermediate steps of the calculation. In the case of 
less inclusive jet observables, one cannot take advantage of the cancel- 
lation of soft and collinear divergences at the integrand level and, at 
present, no systematic method is available to handle these divergences 
at NNLO. Even once the necessary two-loop matrix elements for sev- 
eral processes are calculated, the amount of work needed to provide 
a numerical implementation will be enormous. The main features of 
the dipole formalism, which permit a universal treatment of soft and 
collinear singularities at NLO, seem particularly suited to set up a 
general method for carrying out NNLO QCD calculations. 
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